The Casimir Effect for Generalized Piston Geometries 
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In this paper we study the Casimir energy and force for generalized pistons constructed from 
warped product manifolds of the type / Xy A^ where / = [a, b] is an interval of the real line and A^ 
is a smooth compact Riemannian manifold either with or without boundary. The piston geometry 
is obtained by dividing the warped product manifold into two regions separated by the cross section 
positioned at R e {a, b). By exploiting zeta function regularization techniques we provide formulas 
for the Casimir energy and force involving the arbitrary warping function / and base manifold A^. 

I. INTRODUCTION 






> 

(N 

en 

o 

(N 



X 



The influence that external conditions like electro-magnetic and gravitational fields, non-trivial topology 
and the presence of boundaries have on the Casimir effect has been a very actively considered subject for 



several decades llll-ll4ll. Whereas most research focuses on the impact one particular influence has, due 



mainly to technical progress the combined effect of several aspects is also being analyzed. In particular, 



flat or curved boundaries placed in a curved space-time have received considerable attention recently |115 - 



ISh . Examples are piston configurations, lll9l - l23ll which have the advantage of rendering the Casimir force 



2611 ■ Due to a different 




finite, 112411 although in a curved space-time this is not a generic feature anymore 11251 
geometry in both chambers a complete cancelation of singularities in the force does not necessarily occur 
and ambiguities might remain. However, under certain conditions unambiguous answers can stiU be found 
as was shown recently for conical pistons 

The current article is a continuation of II25L 12611 . where the conical manifold is replaced by a more general 
warped manifold. To analyze the Casimir energy and force we will use the zeta function regularizationll27l. 
2811 and in Section 2 we describe the relevant spectral problem and the geometry of the piston. In Section 
3 we find the analytical continuation of the associated zeta function; part of the construction involves the 
derivation of uniform asymptotics of solutions of initial value problems of an ordinary differential equation 
using the WKB approximation I29l .l30[l . The analytical continuation of the zeta function is the basis for the 
Casimir energy and force results in Section 4. Although our focus is the force on the piston, we present 
results for the energy for each chamber from which the force on the bounding plates of the chambers can be 
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derived. Combining both chambers significant simplifications occur. In the Conclusions we summarize our 
results and outline possible further studies along the lines of our article. 

II. THE GENERALIZED PISTON GEOMETRY AND ZETA FUNCTION 

Let I = [a,b] (Z R. and N be a. J-dimensional smooth, compact Riemannian manifold either with or 
without boundary dN. We consider the bounded manifold M of dimension D = d + I constructed as a 
warped product M = I Xf N assuming that the warping function /(r) > for r € / and that / € C°°(/). The 
local geometry of the warped product manifold M is described by the line element 

ds^ = dr^ + f{r)dl? , (2.1) 

where dl? represents the line element on the manifold N. It is possible to construct a piston configuration 
modeled by the warped manifold M as follows (cf. 11251 12611): Let Nr be the cross section of the manifold 



M positioned at the point R € {a,b). The manifold Nr divides M into two separate regions denoted by 
region / and region //. Region / and region // are represented, respectively, by the D-dimensional compact 
manifolds M/ - [a, R\ X/ A'^ and M// - {R, b] Xf N having boundary dMj = Na U Nr and dMu = NrU Nb- 
It is clear from this construction that the warped product manifold M is given by the union of Mj and Mjj 
along their common boundary Nr, namely M = Mj Ua^^ Mjj. The configuration described is a generalized 
piston where the piston itself is modeled by the cross section A ^^. 



28h in order to obtain the Casimir energy 



In this work we use zeta function regularization techniques 11271 1: 
and force for the generalized piston configuration. For this reason we are interested in the analysis of the 
following eigenvalue equation 

- Am<^ ^ a^(f> , (2.2) 

where A^ represents the Laplace operator acting on scalar functions cp € £}{M). On a manifold with metric 
(12. IK Am is expressed as 

/ d^ /'(r) d 1 \ 

where Aa^ is the Laplace operator on the manifold A'^. Since the fields in one region are independent from 
the fields in the other region, the two spectral problems are independent. This implies that we have a set of 
eigenvalues for each region, namely aj for region / and ajj for region // with corresponding eigenfunctions 
(fij and (pjj. In both regions a solution to the eigenvalue equation (12.21 ) can be written as a product cfip = 
UqM^ ^)^i^)^ with p denoting either / or //, where 0(a>) are the harmonics on N satisfying 

- Aa,(D(w) - y^a)(w) , (2.4) 



and UaJr, v) denotes the solution of 



fir) d 



— r + d^^—^-^ — +ai — I Ua (r, v) = 

dr2 fir) dr P p{r)} " 



(2.5) 



The previous remarks suggest that the spectral zeta function associated with the piston configuration 
acquires contributions from the spectral problems in both regions and can, therefore, be written in the form 



Ms) = ^lis) + ^iiis) , 



(2.6) 



where (,p{s), with p - (I, II), is the spectral zeta function corresponding to either region / or // defined as 



^p(^) - Z 



a 



-2s 



(2.7) 



The explicit dependence of ( 12.51 ) on the eigenvalues y^ of -A/v suggests that the spectral problems on the 
warped manifold M and on the manifold A'^ are closely related. For this reason (m{s) will be expressed in 



terms of the spectral zeta function of the cross section A^ 113 iL 



M 



(Nis)-Y^div)v-^\ 



(2.8) 



where d{v) represents the finite degeneracy of each eigenvalue. From the spectral zeta function (I2.61 l the 



Casimir energy is obtained as follows 02 



m 



,2b 



£cas = lim— -^M s-t) , 



(2.9) 



where ju represents an arbitrary mass parameter. Since ^iwi^) will generally present a pole at s = -1/2, 
performing the limit in (|2.91 l leads to the expression 



1 



1\ 1 /I 



Ec, = :^FP^m[-^\ + :^\-+lnij'\Res^M\-:^\ + 0(e) , 



1 



(2.10) 



where Res and FP denote the residue and the fin ite p art. For piston configurations the Casimir energy 



depends explicitly on the position R of the piston 12411 and the Casimir force can be obtained through the 
relation 



^Cas(^) - --^EcafiiR) 

oR 



(2.11) 



It is not difficult to realize, from the expressions (I2.101 i and (12.1 It . that the Casimir force on the piston is 
well defined only when Res ^m(-1/2) does not depend on the position R. 



III. ANALYTIC CONTINUATION OF THE SPECTRAL ZETA FUNCTION 

The spectral zeta function in (I2.6l l is well defined in the region %(s) > D/2 and in order to compute 
the Casimir energy we need to find its analytic continuation to a neighborhood of s = -1 12. Although the 
eigenvalues a,- are not explicitly known for a general warping function /, they can be found implicitly once 
appropriate boundary conditions are imposed. For definiteness we impose Dirichlet boundary conditions in 
both regions / and //, namely 

Ua,{a,v) = Uai{R,v) = , and Ua„{R,v) = Uai,{b,v) = . (3.1) 

Instead of considering the boundary value problem (I2.5l l and (13.11) . it is convenient to analyze the following 



initial value problem 113 3h 



with pp € C, 



^ fir) d . v^ , 



Upi{a,v) = 0, u'p^(a,v) = 1 , (3.3) 



in region /, and 

Mp„(/?,v)-0, «;„(/?, v)-l, (3.4) 

in region //. The two sets of eigenvalues ai and an are then obtained implicitly as solutions of the equations 

Up,(R,v) = 0, up,,{b,v) = 0. (3.5) 

To avoid displaying similar explicit expressions valid in region / and // we will employ from now on the 

notation Xp with the meaning xj = R and xjj = b. Bearing in mind the last remark and thanks to the relations 

3.11) the spectral zeta functions for region / and // can be rewritten in terms of the integral representation 



1 /^ r\ 

where y^ is a contour that encircles, in the counterclockwise direction, the real zeroes of Uk{xp,v). The 
spectral parameter m has been introduced for technical reasons and it will be sent to zero in the final results. 
By first deforming the contour jp to the imaginary axis and by then performing the change of variable 
k = Z.V, the expression (I3.6l l can be rewritten as 

^p(s) = Y,d{yK;is) , (3.7) 



where 



sin;r5' r°° 99 1 , d 

C{s) = dziv'z^ -m^)-'-\nuUxp,v) ■ (3-8) 

'^ n Jsi dz 

The integral in (|3.8l l is well defined for 1/2 < %is) < 1. In order to perform the analytic continuation 
to a neighborhood of s = -1/2 we will add and subtract, from (|3.8l l. a suitable number of terms from the 
uniform asymptotic expansion, as v ^ 00 and z = k/v fixed, of the functions M,v^(xp, v). 

The desired uniform asymptotic expansion can be obtained by employing a WKB approximation 13411 . 
The starting point of this method is the use of the following expression for the solution of (|3.2| i 



d r 

"2I ' 



M,v,(r,y)-exp -- [/(?) df ^,(z, r) , (3.9) 



with U{t) = fit)/ fit). Once (|T9l ) is substituted into (1X21 ) we obtain 



^+qiv,Z,r)\^Az,r)^0, (3.10) 



with the introduction of the function 

qiv,z,r) = Viv,z,r)-^U'{r)-^U^{r) 



'T-i 



2, 2 , 1 \ ^/"W ^(^-2)/'» 



In order to proceed with the WKB method, we define the auxiliary function 

5(y,z,r) = -ln1',(z,r), (3.12) 

or 

which can be easily shown to satisfy the non-linear differential equation 

S'(y,z,r) - -q{y,z,r)-S^{y,z,r) . (3.13) 

The uniform asymptotic expansion for the function S{y,z,r) as v -^ 00, and hence for *I'y(z, r) through 
(13.121 ). is obtained by substituting the ansatz 

Siv,z,r) ~ vS.iiz,r) + Soiz,r) + ^ ^^ , (3-14) 

1=1 ^ 

into (13.131 ) and by equating like powers of v. This procedure leads to the following expressions for the 
functions Si{z, r) with / = {-1, 0, 1, • • • ), 



/ 1 1 /^ 

5fi(z,r)-±Jz2 + -— , So{z,r) = ---\nSti{z,r), (3.15) 



St{z,r) = - ^ 



2St^{z,r) 



df"{r) did-2)f\r) . 

-27w — ^7w^^ofe^)-^ofe'') 



(3.16) 



and, for the higher asymptotic orders with / > 1, 

1 



Sti(z,r) = -—^ 



2Stiz,r) 



S'r{z,r) + Y^St;{z,r)St,^{z,r) 



«=o 



(3.17) 



Depending on the choice of sign for the functions S i(z, r) one obtains either the exponentially decaying or 
the exponentially increasing part of the asymptotic expansion ior^{z, r). The correct asymptotic expres- 



sion for Uiy^{r, v) contains a linear combination of both behaviors ll34ll and it reads 



+ Bexpi I S'{v,z,t)dt 



Aexp 



Ja 



t)At 



(3.18) 



where A and B are uniquely determined once the initial conditions are imposed. 

By imposing the initial condition (13.31) in region /, (13.41 ) in region 77, and by disregarding exponentially 
small terms we obtain the expansions 



U{t)dt 



in region 7 and 



d r^ 

\nuiv^{R,v) = -ln[S'^{v,z,a)-S {v,z,a)]-- I 

^ Ja 

+ i S'-{v,z,t)dt , 

Ja 



In Uiy,{b,v) - - In [S^{v,z,R)-S-{v,z,R)]-i; f U{t)dt 

2 Jr 

+ r S^{v,zj)dt , 
Jr 



(3.19) 



(3.20) 



in region 77. 

Let us observe at this point that when adding up these two expressions, most terms will become 7?- 
independent. Therefore no contribution to the force (12. lit on the piston will result. 

For the first term on the left hand side of (|3.19l l one has an expansion of the form 



In [S^(v, z, a) - S-(v, z, a)] = In (2v) + ^ In [z^ + -^ I + 



Z 



mz,a) 



i+l 



(3.21) 



where the D,(z, a) are determined through the relation 



In 



1 + tU' + 



f\a) 



oJkiz, a) 



y OJkKZ, a) _ y 



k=\ 



i)i{z,a) 



,,!+l 



(3.22) 



i=l 



with coiiz, a) = S ^(z, a) - S ■ (z, a). Obviously, the same expansion holds for the term 
In [S'^iv, z, R) - S~{v, z, R)] in (13.201 ) once a is replaced with 7?. 



By utilizing the asymptotic expansion (13.141 ). the explicit expression for U{t), and the result (13.211 ) the 
expansion (13.191 ) can be rewritten as 



In Uiy,(R,v) = -ln(2v)-^ln(z2 + ^^) + lln 
J-1 f{a) 



l+z^fia) 



where 



M 



fid) r^ °° 

i{z,a,R)= f S+(z,Odf-£»;-i(z,a), 

Ja 



\+zHKR) 

Mi{z,a,R) 



V' 



(3.23) 



(3.24) 



for / > 1 and Do(z, a) - 0. The asymptotic expansion for In Uiyy{b, v) in region // is obtained from (13.231 ) 
once the replacement a ^ R and R ^ bi& performed. 

Since the process of analytic continuation of the spectral zeta function involves the computation of 
integrals with respect to the variable z of the terms of the uniform asymptotic expansion of the functions 
Uiy^{xp, v), it is convenient to render manifest the dependence on z of the functions Sl'{z, t) and At,(z, a,R). 
From the expression (|3.141 i and the results (|3.15) - (I3.17) one can prove that the functions 5^(z, have the 
general form 



^,^(Z' r) = Y, FkAr)f'^'ir) (l + z^f\r)) 



2k+i 

2 



(3.25) 



k=0 



From the expression for ^0(2, r) and S'^{z, r) in (13.151) and (13.161) it is not difficult to obtain 

FQ,o{r) = 0, Fi^o{r) = ——-, Fo^iir)^ -—— + -— - , 

2/3 (r) 4 /(r) 8 /2(r) 



(3.26) 



^l,lW = T 



3 f'\r) 1 /"(r) 



F2,i{r) = -z 



5 f'\r) 



(3.27) 



4/V) 4/3(r) ' '''''■' 8/6(r) 

The functions F^jir) of higher order are found by using the general form (I3.25t of S'j'iz, t) into the recur- 
rence relation (13.141 ). More explicitly one has, for / > 1, 



Fkj+iir) = -- 



(2k + / - 2):^m;^ + F', .(r)H(i +l-k) + n, iir) 

J in 



where H(x) denotes the Heaviside function and 'Kkjir) is defined through 

(■ m+l i'-m+l 

Z Z Z Fk,,n{r)Fjj.^{r)f''^^i^\r)(\+z^f\r)) 

m=0 k=0 j=0 

i+2 

-J]%,Ar)f"^'ir)(l+z^f\r)) 
(1=0 



Ik+lj+i 
2 



211+1 
2 



(3.28) 



(3.29) 
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The functions Mi(z,a,R) defined in (13.241 ) contain the temis Di(z,a) which, according to the cumulant 
expansion (I3.22I ). have the form 

2( 



n2i^i(z,a) - J]n,,,(a)/'^''(a)(l +zV'(a)) 



-k-i 



(3.30) 



k=0 



with Q,t i(a) found from the relation 



In 



1 + 2 ^7 Z F,,2j-:(a)f''^'\a) (l + zV'(«))""'' 

7=1 ^ \n=Q 

oo 2i 



(3.31) 



1=1 k=0 

The relation (13.301 ) together with (|3.16l l and (13.171 ) allow us to derive the following expressions for 
Aiiiz, a,R): when / = 2ni + 1 and m € No, we have 

2m+2 f^ji 



Mlm+liZ- 



a,R)= Y, \ 

k=o ^" 



dt Ft, 2m+i (0/'^""^' (0 ( 1 + z'/(0)" 



(3.32) 



while for / - 2m and m € N^ we obtain 



2m+l ^R 



M2,niZ,a,R) = Yj [ ^tFk,2m{t)f^^^'^{t)(\+Z^f\t)) 
k=Q ^" 

2m 

- Y,^t,n,ia)f'^""ia)(\+z'f\a))' 



k-m 



-k-m 



(3.33) 



k=0 



We would like to remind the reader that similar results for A1,(z, R, b) are obtained in region // once the 
replacement a ^ R and R ^ bis made. 

Since the relevant uniform asymptotic expansion is now completely known we can proceed with the 
analytic continuation of the spectral zeta function (I2.6l l. By adding and subtracting in (13.81 ) N = D terms of 
the uniform asymptotic expansion (13.231 1 we obtain, in region /, 

D 



(3.34) 



(i(s,a,R) = Zj(s,a,R) + ^ Af\s,a,R) . 
The function Zj{s,a,R), in the massless case, has the following integral representation valid for %{s) > -1 

Zi{s, a,R) = ^^^^ y d{v)v'^' Azz~^' — \ In m^^iR, v) + In (2y) 
^ V Jo dz{ 



+ xln U^ + 



1 



na)j 4 



-lln 



l+zV» 
l+z^p{R) 



2 " /(/?) 






(3.35) 



The contributions A {s, a, R) are obtained, instead, by integrating with respect to the variable z the asymp- 



totic terms in (13.231 ) and they read, in the massless case. 



iW 



A'_[{s,a,R) 



1 r(.-i) 
2 V^ r{s) 



^n\s 



w 



f'-\t)dt , 



A\;\s,a,R) = -^(^(s)[f%a)+f\R)] , 



(3.36) 



(3.37) 



for n € No, and 



A2,i(^,a,/?) 



r(^) 



I \^ 2n+2 y(s + j + n + 



(N\s + n + -]y 



i) 



;=o 



{j + n + y 



X Fj,2n+l{t)f 



-2(s+j+n)+l 



{t)dt , 



^ 2h+1 p ^ 



r( ^ + 7 + ?i) 



X 






+^+"'(Odf - a;,„(a)r<-^+-'+")(«) 



(3.38) 



(3.39) 



when « € N""" and where Q.2n+i,n - 0. In region // we obtain similar results for the analytic continuation of 
the spectral zeta function. In fact, one has 

D 



Cn{s,R,b) = Zu{s,R,b) + Y^ Af\s,R,b) . 



(3.40) 



(•=-1 



The terms Zjj{s, R, b) and A* {s, R, b) can be obtained from, respectively, (13.351 ) and (I3.361 l- (I3.39I ) once we 
replace a ^ R and R ^ b. 

IV. THE CASIMIR ENERGY AND FORCE FOR GENERALIZED PISTONS 



In order to compute the Casimir energy in the framework of zeta function regularization we employ the 
relation (12.101 ). By setting 5 = e - 1/2 and by taking into account the meromorphic structure of ^n{s) iSl] 
we obtain for the asymptotic terms (I3.361 l- (l3.391 l in region / the expansions 



r\t)lnf{t)dt 



A«(£-l/2,a,/?) = — ^^(-1)J /-2(f)df + _^^(-i)J ^ 

1 r^ 

+ — [4(-l) + (ln4-l)^^(-l)]j r\t)dt^O{s), 



(4.1) 



10 



and 



Af{s-ll2,a,R) = -lRes^;v(-^j[rHa) + rH^)] 
-^Res^;v (-^) [r '(a)ln/(a) + rk/?)ln/(/?)] + 0(£) , 



+ ^^n{0) 



-R 



R 



if Fo,i(01n/(Od? + (ln2-l) | Fo,i(Od? 

2 ^. .^ ^R 






The remaining asymptotic terms with n > \ have the form 

^ 2«+2 



W + n) 



A2^,is-l/2,a,R) = — Res^^(«) V —^ ^— 



:nj,n{a,R) 



+ Sj,n(a,R) 



and 



2n 



1 / i\2^ir(j + «-i) 

(e- l/2,a,/?) - — Res^^ « - - V ^„^ . '^ 



C;,„(a,/?) 



+r;>(a,/?) 



i=0 



where for typographical convenience we have introduced the functions 



1 



:Hj,„{a,R) = - I Fj^2n+i(t)f^^^"mt , 



(4.2) 



(4.3) 



1 2n+2 J-,, . , 

+ — FP ^N (n) V , ^^^"\. Jlj,n{a, R) + 0(e) , (4.4) 



(4.5) 



(4.6) 



pR 

3j,n{a,R) - Fj,2n^i{t)f^^^\t)lnf{t)dt + (Hn^j^i - 2 + 2\n2) J{j^n{a,R) , (4.7) 



1 r^ 

Cj,n(a, R) = 2 ^i. 2„(0/^^^"~k0dr - Tj,n(a) , 



with 



(4.8) 



Tj^nia) = -Qj,n{a)f'^^"'\a) , 



(4.9) 



11 



and 



Tj,n{a, R) = Fj^ 2„(0/'^'""k0 In mdt - Qy, nia)f'^^''-\a) In /(a) 

Ja 



+ 



n+j-1 

Z 

k=0 



Ik -I 



Cj,n{a,R), 



(4.10) 



where H„ represent the harmonic numbers. In addition, since Zj{s,a,R) is an analytic function in the 
halfplane %is) > -1 the value s = -1/2 can simply be set in (13.351) . Once again, the same results are 
obtained in region II with the replacement indicated before. 

From the results (I4.1I )- (I4.5I ) it is not very difficult to compute the Casimir energy associated with region 
/. In fact, by utilizing the formula (12.101 ). with the spectral zeta function ^m replaced by ^/, we obtain the 
following expression 



2 \ 2 



E'^^Ja,R) = ^Z,{-^,a,R\ + -Q{a,R) 



V- 1 2n+2 



2V^ ^ ^ r(j + n + ^) 
+ Res ^N in) Sj,n{a,R)] 



\j + n) 



+ Res^N\n- -\Tj,n{a,R) 



VV^N\n--\Cj,n{a,R) 



^\[\^W 



^^;v(-i) J r\t)&t 

+ ^^yv(0) J^ Fo,i(Odf - ^Res ^n (-^) [r\a) + f-\R)\ 
+ -p y Res^w(n)y — ^ 



m+n) 



+ n + 



D 



:Rj,n{a,R) 



(4.11) 



1 



[f] 



^^2^T(i + n-\) 



+ ^y Res^/v L-- y — ^ ^-^Cin{a,R)\ + 0{E) 

where the notation [.i] stands for the integer part of x and, for convenience, we have defined the function 



1 r^ 

Q{a,R) = — [^(-1) + (ln4 - l)^;v(-l)] J /"'(Odf 



+ 2;^^w(-l) I ri01n/(0d? 



^Res^w(-i][/-i(a)ln/(a) + /-!(/?) In /(/?)] 



12 



+ :^^n{0) 



Ja 



Fo,i{t)\nf{t)At + {\n2-\) Fo,iit)dt 



f 

Ja 



:^^N{0)i]-f^ f Fj,,{t)fJ{t)dt 



(4.12) 



The Casimir energy in region // can be obtained from (14.111 ) once the replacement a ^ R and R ^ b is 
performed and by considering Zjj{s,R, b) in place of Z/(i^, a,R). 

In order to compute the Casimir force on the piston we need, according to (12.101 ) and (12. lib , the residue 
and finite part of the spectral zeta function Cm{s) at 5 = -1/2. By adding the contributions coming from 
region / and region // we obtain, for the residue, the following expression 

- lRes^A,|-iJ [f-\a) + f-\b) + 2/-i(/?)] 
1 fS ^^ Yij + n) 



[?] 



+ > Res/ft, n-- > — ^ ' 

V^^ I 2JZ^ r( 

' (1=1 ^ ' 7=0 

X [Cj,n{a,b)-Tj,n{R)] . 



\i + n) 



Furthermore, for the finite part we have 



FP^mI-^) = Zii-^,a,R\ + Zn(-\,R,b\ + Q{a,b) + }f{R) 



. 2 2;i+2 -re • , \ 



+ ^ y y -r^ '—[¥P(N{n):Rj,n{a,b) 

+ Res^N{n)Sj^nia,b)] 

„=i,=o ro + n) [ \ 

+ Res ^N In - ^ J [rj,n{a, b) - Xj,n{R)) 
where we have introduced the functions 



1 [f]2«.ir(j + „-i)r ^ 



}f{R) = -^FP^N(-jjr\R)-Res^N(-\]r\R)inf{R) , 



and 



Xj^R) = aj,niR)f'^^"-\R) 



n+j-l 



In 



*:=0 



(4.13) 



(4.14) 



(4.15) 



(4.16) 
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Thanks to the formulas (14.131 1 and (14.141 ) the Casimir force on the piston, which is the quantity of main 
interest, can be evaluated according to (12. lit . By differentiating the above results with respect to the position 
R of the piston and by denoting with a prime such derivative we obtain 

FcsiR) = -^Z'J-^,a,R]-jZ'jJ-j,R,b\-^}/'{R) 



with 



1 [|]2^i r(j + n-i) 



FP^M{n--\T;,^(R) 



+ Res^Nin-^\?Cj,{R) 



1 /I 



-ilr^"^1i^-^n-2)m 



1\/'W 



[#] 



n-ir(;>.-i) 



-J^Res^..--!;^^— ^77,„(/.) .0(.), (4.17) 



J/'(/?) 



n=l ^ ' ;=0 



|;||[^FP^.(-l).Res^.(-il(l-ln/W) 



(4.18) 



T-;„(/?) = ^/^■^'"-'(/?) [q;;„(/?)/(/?) + (2« + 2; - l)Qy,„(/?)/'(/?)] , 



(4.19) 



and 



'^^R) = 277,„(/?) 



(7+7-1 



In /(/?)+ 2 



1 



*:=0 



2/t- 1 



+ 2Tj,„{R) 



f'JR) 
fiR) 



(4.20) 



We would like to make a few remarks regarding the result (14.171 ) for the Casimir force on the piston. The 
term proportional to (1/e + In/j.^) represents the ambiguity which is generally present in the Casimir force. 
The ambiguity is proportional to the residue of ^/v(i') at the points s = n - 1/2 with < n < [D/2]. The 
residues are related to the heat kernel coefficients ao/i-n associated with the manifold N which shows that 
the ambiguity is purely of geometric nature. It is clear, from the previous remark, that the force on the 
piston becomes a well defined quantity when the manifold N is even-dimensional and without boundary, a 
result that also holds for Casimir pistons modeled by a generalized cone 125112611 . We would like to point 
out that the term in expression (14. 17l i responsible for the ambiguity of the force is proportional to the first 
and higher derivatives of the warping functions /(/?). This can be easily seen for the term proportional to 
Res ^N (-1/2). Moreover, for the remaining terms we notice that T'-^{R) contains first and higher derivatives 
of /(/?) through Q.jn(R) and Q.'.{R) (see, for instance, Eqs. (13.311 1 and (I3.28l l). The last remark shows that 
when the warping function is constant, namely the warped manifold M becomes a generalized cylinder, the 
force is always unambiguous regardless of the geometry of the manifold N, as has been observed earlier 
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V. CONCLUSIONS 

In this work we have studied the Casimir energy and force for massless scalar fields endowed with 
Dirichlet boundary conditions in the framework of a generalized piston constructed from a warped product 
manifold. In order to compute the Casimir energy and force in this setting we have performed the analytic 
continuation to a neighborhood of s = -1/2 of the spectral zeta function associated with the Laplace 
operator defined on warped product manifolds. The method employed for the analytic continuation is based 
on the contour integral representation of the zeta function and the WKB expansion of the eigenfunctions of 
the Laplacian on the manifold M. This procedure allows one to study the spectral zeta function associated 
with manifolds posses sing quite general geometries. A more complete and detailed account of the method 



has been presented in B6ll . 

The results that we have obtained for the Casimir energy and force are valid in any dimension D and 
for any smooth, compact piston N. In addition, both the expressions for the energy and force depend 
explicitely on the spectral zeta function ^iv(s) which shows how these quantities are related to the geometry 
and the topology of the piston. The interest in studying generalized piston configurations lies in the fact 
that, for arbitrary warping functions, the two chambers forming the piston do not have the same geometry. 
This implies that the behavior of the Casimir force for generalized piston configurations is different from 
the usual Casimir pistons considered in the literature where the two chambers have the same geometry. 
We would like to mention that the formulas for the Casimir force on the piston, although very general, 
are somewhat implicit. However, once a specific warping function /(r) and manifold N are chosen more 
explicit results can be obtained. 

In fact, this would be an interesting continuation of the presented research. Namely, study in more detail 
particular examples of generalized pistons which are also of relevance in physical situations. Furthermore 
other boundary conditions can be envisaged by essentially the same methods. 
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